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Abstract. We develop a numerical chemical model al- 
lowing for radial flows of gas, with the aim to analyse 
the possible role of gas flows in the chemical evolution of 
the Galactic Disc. The dynamical effects of the Galactic 
Bar on the radial gas profile of the Disc are especially 
addressed. 
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dance gradients - Galaxy: gas distribution - Galaxy: Bar 



1. Introduction 

The observed chemical and spectro-photomctric proper- 
ties of galaxies are one of the main sources of information 
for our understanding of galaxy formation and evolution. 
The corresponding theoretical modelling involves star for- 
mation (SF) as a basic ingredient. Unfortunately, this pro- 
cess is rather poorly known on the large scales relevant 
to galaxy evolution. Portinari & Chiosi (1999, hereinafter 
PC99) analysed the effects of adopting different SF laws 
in a chemical model for the Galactic Disc, a system which 
we can study in great detail. In this paper we address an- 
other phenomenon which can bear interesting effects on 
the chemical evolution of galaxies: radial gas flows. A few 
papers in literature (Section ^ demonstrate that radial 
gas flows influence chemical models for the Disc, espe- 
cially in their predictions on the metallicity gradient. It 
is therefore interesting to discuss the radial profile of the 
Disc with models including also radial flows, in addition to 
various options for the SF law. In particular, radial flows 
can help to overcome some difficulties that "static" mod- 
els find in reproducing, at the same time, the metallicity 
gradient and the radial gas profile of the Disc (PC99). 

We develop a chemical model with radial gas flows as 
a multi-dimensional extension of the model of Portinari 
et al. (1998, hereinafter PCB98). The model is described 
in Section^ and in the appendices. In Section^ we discuss 
the general qualitative effects of superposing radial flows 
upon a chemical model. In Section ^|we present models for 
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the Galactic Disc with radial gas flows and different SF 
laws, showing that radial flows provide an alternative or 
additional dynamical effect to the "inside-out" formation 
scenario to explain the metallicity gradient. Section ^| is 
dedicated to qualitative simulations of the dynamical ef- 
fects of the Galactic Bar upon the gas distribution, with 
the aim to reproduce the molecular ring around 4 kpc, 
which static models cannot account for (PC99). Section [?] 
contains a final summary and conclusions. 



2. Radial flows: previous literature 

The possibility that radial flows play a role in establishing 
the radial metallicity gradients in galactic discs was first 
suggested by Tinsley & Larson (1978). Following Lacey & 
Fall (1985), we mention that radial gas flows in a disc can 
be driven by three main mechanisms: 

1. the infalling gas has a lower angular momentum than 
the circular motions in the disc, and mixing with the 
gas in the disc induces a net radial inflow with a ve- 
locity up to a few km sec -1 ; 

2. viscosity in the gas layer induces radial inflows in the 
inner parts of the disc and outflows in the outer parts, 
with velocities of ~0.1 km sec -1 ; 

3. gravitational interactions between gas and spiral den- 
sity waves lead to large-scale shocks, dissipation and 
therefore radial inflows of gas (or outflows in the outer 
parts) with typical velocities of ~0.3 km sec -1 (e.g. 
Bertin & Lin 1996 and references therein); much larger 
velocities can be achieved in the inner few kpc in the 
presence of a barred potential. 

In summary, radial flows are plausible with velocities of 
~0. 1— 1 km sec -1 , and they are expected to be inflows 
over most of the disc. Observational upper limits permit 
radial inflows in the Galactic Disc with velocities up to 
5 km sec -1 at the present time. For further details, see 
Lacey & Fall (1985) and references therein. 

The first of the above mentioned mechanisms was mod- 
elled in detail by Mayor & Vigroux (1980), and later 
by Pitts & Tayler (1989, 1996), Chamcham & Tayler 
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(1994). The effects of a generic inflow velocity pro- 
file on chemical evolution models has been explored by 
Lacey & Fall (1985), Tosi(1988), Gotz & Koppen (1992), 
Koppen (1994), Edmunds & Greenhow (1995). A dif- 
ferent approach is that of viscous disc models which, 
rather than imposing arbitrary radial velocity patterns, 
describe the evolution of the gas distribution in the disc 
self-consistently, following the model suggested by Lin & 
Pringle (1987). Viscous chemical models have been devel- 
oped by Clarke (1989), Yoshii & Sommer-Larsen (1989) 
and Sommer-Larsen & Yoshii (1990), Thon & Meusinger 
(1998). All these studies show how radial inflows can 
steepen the metallicity gradients with respect to static 
models, especially if an outer cut-off of SF is assumed. 

3. Modelling radial flows 

We formulate our chemical model with radial flows as a 
multi-dimensional extension of the static model of PCB98 
and PC99, an open model where the disc forms gradually 
by accretion of protogalactic gas. The disc is divided in 
N concentric rings or shells; in each ring fc the gaseous 
component and its chemical abundances evolve due to: 

1. depletion by SF, which locks up gas into stars; 

2. stellar ejecta which shed back enriched material to the 
interstellar medium (ISM); 

3. infall of primordial protogalactic gas; 

4. gas exchange with the neighbouring rings because of 
radial flows. 

The set of equations driving the chemical evolution of the 
fc-th shell is: 

f t Gi{r k ,t) = -X l (r k ,t)^{r k ,t) + 



+ J^;^(r k ,t-T M )R l (M)^(M)dM + 



(1) 



inf 



+ [fMn,t)} rf 

where the various symbols are defined here below. 

Primordial gas is accreted at an exponentially decreas- 
ing rate with time-scale r: 



<Jinf{rk,t) = A(r k ) e < r h) 



(2) 



A(r k ) is obtained by imposing that the integrated contri- 
bution of infall up to the present Galactic age to =15 Gyr, 
corresponds to an assumed exponential profile <JA(r k ): 



A(r k ) = 



a A {r Q )e 



(3) 



r(r fe )(l - e-^M^-)) r(r fc )(l - e-^Mn,)) 

Indicating with a g (r k , t) the surface gas density, we define 
the gas fraction: 



G(r k ,t) 



<?A(r k ) 



(4) 




Fig. 1. 

shell 



Scheme of the gas flow through the fc-th model 



and the normalized surface gas density for each chemical 
species i: 



G i (r k ,t)=X i (r k ,t)G(r k ,t) 



(5) 



where Xi is the fractionary abundance by mass of i. 

The 1 st term on the right-hand side of Eq. (Q) repre- 
sents the depletion of species i from the ISM due to star 
formation; see PC99 for the various options concerning the 
SF rate, 4'(r, t). The 2 nd term is the amount of species i 
ejected back to the ISM by dying stars; the returned frac- 
tions Ri{M) are calculated on the base of the detailed 
stellar yields from PCB98 and keep track of finite stellar 
lifetimes (no instantaneous recycling approximation IRA). 
The 3 rd term is the contribution of infall, while the A th 
term describes the effect of radial flows. Full details on 
the first three terms can be found in the original static 
model by PCB98 and PC99. The novelty in Eq. (Q) is the 
radial flow term, which we develop here below. We will 
adopt the simplified notation a gk = a g (r k , t) and the like. 

Let the fc-th shell be defined by the galactocentric ra- 
dius r k , its inner and outer edge being labelled as r k _ i and 
r fe+i- Through these edges, gas flows with velocity v k _i 
and v k+ i, respectively (Fig. |l|). Flow velocities are taken 
positive outward; the case of inflow is correspondingly de- 
scribed by negative velocities. Radial flows through the 
borders, with a flux F(r), contribute to alter the gas sur- 
face density in the fc-th shell according to: 



da, 



dt 



1 



y>2 y->2 



rf 7T C 

The gas flux at r k+ i can be written: 



F(r k+h )~F(r k _ h ) (6) 



F(r k+ i) = 2Trr k+ iv k+ i xK+ \ ) °gk + 1 ) cr 9 (fc+i) (7) 

where x( x ) is the step function: x( x ) = 1 or for 
x > or < 0, respectively. Eq. (0) is a sort of "upwind ap- 
proximation" for the advection term to be included in the 
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model equations (e.g. Press et al. 1986), describing either As to the innermost edge, we assume that the first shell 
inflow or outflow depending on the sign of v k+ i. An anal- is symmetric with respect to r±: 



3ri - r 2 



ogous expression holds for F(r k _ 1 ). 

Let's take the inner edge r k _i at the midpoint be- r 
tween r k _i and r k , and similarly for r k+ i (Fig. hi). Writ- 
ing Eq. (||) separately for each chemical species i, in terms and that vi < always, since we cannot account for out- 
flows from still inner shells, not included in the model. For 
fc = 1, Eq. (H) then becomes: 



of the Gi's we obtain the radial flow term of Eq. ([!]) as 



[■&Gi(r k ,t)] rf = a fc Gi(r fe _!,t) - k Gi(r k ,t) + 



(8) 



where: 



(10) 



rf 



r k + 



~k-l+ r k-\ 



x(f t -')% 



r k -i+r k 
2 r k+1 -r k _ 1 



a A(k-l) 
&Ak 



0k 



Ik 



r k + - 



x(-v k -i)v k _i 



r k -i+r k 
1 ''fc + l+ffc-l 



X(v, 



k+ 



i )v 



r k +r k+ i 



(9) 



r k-l+ r k4 



x(-v k+ i)v k+ 



r k +r k + i 



2 r fc+ i-r fc _i 



&A(k + l) 

&Ak 



The terms on the right-hand side of Eq. (g) evidence the 
contribution of the 3 contiguous shells involved: the first 
term represents the gas being gained in shell k from k— 1, 
the second term is the gas being lost from k to fc-1 and 
fc+1, and the third term is the gas being gained in k from 
fo+1. The coefficients (^) are all > and depend only on 
the shell k, not on the chemical species i considered. If 
the velocity pattern is constant in time, a k , 0k and 7& are 
also constant in time. 

Notice that in the case of static models the final surface 
mass density is completely determined by the assumed ac- 
cretion profile, namely a(r k ,tc) = ita (?"&)■ Therefore, in 
static models the radial profile for accretion can be directly 
chosen so as to match the observed present-day surface 
density in the Disc (see PCB98 and PC99). The inclusion 
of the term of radial gas flows alters the expected final den- 
sity profile and a(r k ,tc) 7^ o-^(r k ). Hence, a(r,tc) cannot 
be assumed in advance and is only known a posteriori (see 
§P ; at the end of each simulation we need to check how 
much radial flows have altered the actual density profile 
c(?"fc; to) with respect to the pure accretion profile <JA(rk)- 
With the slow flow speeds considered (v < 1 km sec -1 ), 
the two profiles will not be too dissimilar anyways. 

3.1. Boundary conditions 

Eq. (||) needs to be slightly modified in the case of the 
innermost and the outermost shell, since the shell k-1 or 
fc+1 are not defined in these two respective cases. 



3.1.1. The innermost shell 

Our disc models will extend down to where the Bulge be- 
comes the dominating Galactic component (ri = 2.5 kpc). 



with: 
0i = - 

I 

71 = - 



1 



3r x - r 2 



r 2 - n 



X(V3)V3 



ri + r 2 
r 2 - ri 



, \ 1 n + r% A2 

-X\— V3 ) V3 



(ii) 



3.1.2. Boundary conditions at the disc edge 



As to the outermost shell (k = N), we need a boundary 
condition for the gas inflowing from the outer disc. We 
assume a SF cut-off in the outer disc, while the gaseous 
layer extends much further. In fact, in external spirals HI 
discs are observed to extend much beyond the optical disc, 
out to 2 or even 3 optical radii. A threshold preventing 
SF beyond a certain radius is expected from gravitational 
stability in fluid discs (Toomre 1964, Quirk 1972) and has 
observational support as well (Kennicutt 1989). For the 
Galactic Disc we assume no SF beyond the last shell at 
rjy = 20 kpc, which is the empirical limit for the opti- 
cal disc and for HII regions and bright blue stars tracing 
active SF. Gas can though flow in from the outer disc; 
extended gas discs might actually provide a much larger 
gas reservoir for star-forming spirals than vertical infall, 
at least at the present time when the gravitational settling 
of the protogalactic cloud is basically over. 

If no SF occurs in the outer disc, the evolution of the 
gas (and total) surface density can be expressed as: 



do- , . . , % _ 
— (r,t) =A(r)e 



r(r] 



Id 

r or 



Vr > r 



(12) 



(e.g. Lacey & Fall 1985, their model equation with the SF 
term dropped). Here, with no SF, a = a g and abundances 
always remain the primordial ones (X iin f). Let's assume 
the following simplifying conditions for the outer disc: 



1. the infall time-scale is uniform: 



r(r) = r(r N ) 



Vr > r 



AM 



2. the inflow velocity is uniform and constant: 



v(r,t) = v N+i 



Vr>r N+ i, Vi 
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the infall profile is flat: 
A{r) = A ext 



3.2. The numerical solution 



Vr > r 



NA 



in accordance with observed extended gas discs in spi- 
rals, showing a much longer scale-length than the stel- 
lar component. 

With these assumptions, Eq. (|l2|) becomes: 

da da . t v ,„„. 
— + v— = Ae ~ - -a (13) 
at or r 

where we indicate r = t(tn), v = v N+ i and A = A ext 
to alleviate the notation. Eq. (|l^) has a straightforward 
analytical solution (Appendix B): 

a(r, t) = At x 

"(l - e -*) + £ (r (e-^ - e~$ ) - (t - T rf ) 



(14) 



where T r f > is the time when radial inflows are as- 
sumed to activate. Eq. (Q) is our boundary condition at 
the outermost edge. 

Notice that dl4| ) is the solution of (JL3|) in the ideal- 
ized case of an infinite, flat gas layer extending boundless 
to any r > r^ (see also Appendix B). Of course, this 
does not correspond to gaseous discs surrounding real spi- 
rals; but since we will consider only slow inflow velocities 
(v < 1 km sec -1 ), with typical values of rjv = 20 kpc and 
to — 15 Gyr, the gas actually drifting into the model 
disc shells will be just the gas originally accreted within 
r ~ 35 kpc. Therefore, the boundary condition ( |l4| ) re- 
mains valid as long as the gas layer stretches out to 
~ 35 kpc, a very plausible assumption since observed 
gaseous discs extend over a few tens or even ~ 100 kpc. 

3.1.3. The outermost shell 

We take a reference external radius r ext > r/v in the outer 
disc where the (total and gas) surface density a{r ext ,t) = 
&ext(t) is given by the boundary condition (|l4| ); typically, 
r N — 20 kpc and r ext ~ 21 kpc. We take the outer edge 
of the shell at the midpoint: 



rN + r ex t 



and replace 

Ai(fe+i) a g (k+i) — y Xij n f a ex t 

in Eqs. (||) and (0), since the primordial abundances Xi^ n f 
remain unaltered in the outer disc, in the absence of SF. 
We thus write the radial flow term for the iV-th shell as: 



[ftGi{TN,t)] rf 



&n Gi(r N -i,t) 

+ UJ t (t) 



f3 N Gi(r N ,t) + 



(15) 



where: 



-X Unf x(-«i\r+i) v N+i) x 



(16) 



Using (g), (10) and (p^[), the basic set of equations (Q) can 
be written as: 



0i (t) Gi (n , t) + 7l Gi (r 2 , t) + Wi (n , t) 



f t G t (r k ,t) = a k Gi(r k -i,t)+# k (t)Gi(r k ,t)+ 
+j k Gi(r k ,t) + Wi(r k ,t) 

k = 2, ...TV - 1 

f t Gi(r N ,t) = a N G l (r N ^ 1 ,t) + d N (t)G l (r N ,t)+ 
+Wi(r N ,t)+LJi(t) 

where we have introduced: 

fc (t) = - (r)(r k ,t) + f3 k ) <0 

r](r k ,t) ee §(r fc ,t) 

Wi(r k ,t) ee f™;V(r k ,t-T M )Ri(M)${M)dM + 



inf 



We refer to PCB98 for further details on the quantities rj 
and Wi, appearing also in the original static model. Ne- 
glecting, for the time being, that the rj's and the Wi's 
contain the GVs themselves, we are dealing with a linear, 
first order, non homogeneous system of differential equa- 
tions with non constant coefficients, of the kind: 



dt 



= A(t)Gi(t) + Wi(t) 



(17) 



There is a system (|T^) for each chemical species i , but the 
matrix of the coefficients A(t) is independent of i. 

We solve the system by the same numerical method 
used for the original equation of the static model - 
see Talbot & Arnett (1971) and PCB98 for details. We 
just need to extend the method to the present multi- 
dimensional case (^) . If we consider the evolution of the 
GVs over a short enough timestep t% — to = At, the various 
quantities rj(r k ,t), i?fe(t) and Wi(r k , t) will remain roughly 
constant within At; similarly to the method for the static 
model (see PCB98), within At we approximate them with 
the values rj k , 0^ and Wi{r k ) they assume at the midstep 
tx = to + hAt. Over At, (Jl7|) can then be considered a 



system with constant coefficients A = A(ti), and Wi(t) 
becomes a constant vector, which allows for the analytical 
solution: 



G i (t 1 ) = e AtA G i (t ) 



to 



where e*" 4 indicates the matrix: 



(e Al *«i 



e XNt u N ) («x 



Wi 



U N , 



(18) 



(19) 



r N -i+2r N +r ext r £xt -r N - 



with Afc the eigenvalues of A and u k the corresponding 
eigenvectors. The matrix e tA and the explicit expression 
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of the solution (|l^) for the £?,'s are calculated in Ap- 
pendix A, resulting in: 



Gi(r k ,ti) = a k x 

x 



-(5>*-l+jS h _i)At 



Gi(r fc _i,< ) + 



-7fc 



P) fe -i+/3fc-i)-p? fe +A=) 

+ e-<"*+^ A * Gi(r fc ,t ) + 

(~fc + l+/3* ! + l)At_ e -(- fc +^ fc )At 



-(lJfc-l+/3 fc _i)At 



Vk+Pk 

Wi{r k ) i 



'jfc-i+^fc- 

(~fc+<3fc)At 



(20) 



+ 



+ 7fe 



(v k +£k)-(vk+i+0k + i) 

.g-tlfc + l+^fe+lJAt 1 



X 



"(^fc+ffcJAt 



n k+1 +Pk+i 



where we intend a\ = 0, 7^ 
shell one should replace: 



■ik+Pk 

0, and for the outermost 



Wi(r N ) — ► Wi{r N ) +uji = W l {r N ) + 



1 

At 



uii(t) dt 



uii being defined by (|T6|). 

Similarly to the case of an isolated shell (see PCB98) 
the system ( p0| ) docs not provide the final solution for 
Gi(r k ,ti), since the fj^s and the Wi's on the right-hand 
side actually depend on the G»'s; it just represents a set 
of implicit non-linear expressions in Gj(rfe,ti). As in the 
original static model, we can neglect the dependence of 
Wi{r k ) on the Gi's and consider only that of rj k (see 
PCB98 for a detailed discussion). We must finally find the 
roots of the system (20) by applying the Newton- Raphson 
method, generalized to many dimensions (cfr. Press et al., 
1986). Such a system holds for each of the chemical species 
i considered, so at each iteration we actually need to solve 
as many systems as the species included in the model. Full 
details on the mathematical development of the model and 
its numerical solution can be found in the appendices and 
in Portinari (1998). 

We tested the code against suitable analytical coun- 
terparts and obtained the following conditions for model 
consistency (Appendix B). 

1. Rather small timesteps are needed for the numeri- 
cal model to keep stable; the required timesteps get 
smaller and smaller the higher the flow velocities con- 
sidered, and the thinner the shells. 

2. To describe gas flows in a disc with an exponential den- 
sity profile, the shells should be equispaced in the loga- 
rithmic, rather than linear, scale (so that they roughly 
have the same mass, rather than the same width). 

We modelled the Galactic Disc using 35 shells from 2.5 
to 20 kpc, equally spaced in the logarithmic scale, their 
width ranging from ~ 0.2 kpc for the inner shells to ~ 
1 kpc for the outermost ones. With such a grid spacing, 
and velocities up to ~ 1 km sec , suitable timesteps are 
of 10~ 4 Gyr (Appendix B; see also Thon & Meusinger 



1998). This means that roughly 1.5 x 10 5 timesteps, times 
35 shells, are needed to complete each model, which would 
translate in excessive computational times. This drawback 
was avoided by separating the time-scales in the code. 

1. The timestep At used to update the "chemical" vari- 
ables (77, Wi, etc.) is the minimum among: Ati which 
guarantees that the relative variations of the Gt 's are 
lower than a fixed e; At2 which guarantees that the 
total surface mass density a{t) increases by no more 
than 5%; At3 which is twice the previous timestep of 
the model, to speed up the computation when pos- 
sible; At4 which guarantees the Courant condition 
At < v Ar, indispensable for the stability of a numer- 
ical algorithm describing flows. So, At is basically set 
by the requirement that the chemical quantities do not 
vary too much within it, and it can get relatively large 
(up to 0.2 Gyr), especially at late ages when the vari- 
ous chemical variables evolve slowly. 

2. It is only the numerical solution ( |2C| ) which needs 
very short timesteps to keep stable. Therefore, 
once the chemical variables are upgraded, the main 
timestep At is subdivided in much shorter timesteps 
St = 10~ 4 Gyr, upon which the solution ( |2C| ) and its 
Newton-Raphson iteration are successively applied to 
cover the whole At. Only then a new upgrade of all 
the rj k 's and Wi's is performed. 

This trick keeps the code roughly as fast as if it would 
evolve with a single time-scale At, and yet it gives the 
same results as the "slow" version in which all quantities 
are upgraded at each St = 10 -4 Gyr. 

Table 1. Parameter values and resulting metallicity gra- 
dients for models S15RF, O10RF and DRRF 



model 



v(r) 



d[Q/H] 
dr 



S15a 


4 


0.35 


0.2 


3 





-0.03 


S15RFa 


4 


0.35 


0.2 


3 


-1 


-0.053 


S15RFb 


5 


0.4 


0.32 


3 


-1 


-0.047 


S15RFc 


4 


0.35 


0.2 


3 


\ V N+i = 


-0.083 


S15RFd 


7 


0.35 


0.35 


3 


/ ^ 

1 V N+l = 


-0.073 


S15RFe 


6 


0.37 


0.25 


3 


Fig. ^ 


-0.063 



OlOa 


4 


0.19 


0.2 


3 





-0.03 


OlORFe 


7 


0.2 


0.25 


3 


Fig. \ 


-0.07 




DRa 


4 


0.42 


0.2 


3 





-0.07 (r > r ) 
flat (r < r ) 


DRRFe 


4 


0.5 


0.27 


3 


Fig. ? 


-0.059 
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Fig. 2. Left panels: final profiles of the total surface density from models S15RF with inflow from the outer disc, in 
linear (upper panel) or logarithmic (lower panel) plot; the reference observed profile of models with scale length 4 kpc, 
as adopted in models with no flows (PC99), is shown as a thick solid line. Upper right panel: radial gas density profiles 
compared to the observed one. Lower right panel: predicted radial metallicity profiles compared to the observational 
data (data and symbols as in PC99). 



4. Effects of radial inflows on the disc 

After presenting our new model with radial flows, we in- 
vestigate their general effects on the chemical evolution of 
the disc. Here we first analyse how the predicted metallic- 
ity gradient and gas and total surface density profiles of 
a static model are altered by a superimposed uniform gas 
inflow, while in Section ^ we will suitably combine radial 
flows with various SF laws to match the observed radial 
profile of the Disc. The characteristics and parameters of 
the various models presented in this and in the next sec- 
tion are summarized in Table [y. 

To gain a qualitative understanding of the effects of 
radial gas inflows on chemical evolution, we impose a uni- 
form inflow pattern v — —1 km sec -1 upon a static chemi- 
cal model and compare the new outcome with the original 
static case. Radial flows over the disc are expected to be 
mainly inflows, with velocities from 0.1 to ~1 km sec -1 
(see §||); imposing a uniform inflow of 1 km sec -1 is there- 



fore a sort of "extreme case" , considered here for the sake 
of qualitative analysis. Anyways, previous studies in liter- 
ature suggest that the effects of radial flows saturate for 
much higher velocities (Koppen 1994). We will consider 
both the case of inflow from the outer gaseous disc and 
not (v N+ i = —1 or 0, respectively). 

All models with radial flows are rescaled so that the 
final surface density at the Solar ring (r© = 8.5 kpc) 
corresponds to 50 M© pc -2 . Namely, with radial flows 
c{tq, to) 7^ <ja(?o) and it cannot be imposed as an input 
datum (see §^), but the zero-point A(r©) of the exponen- 
tial accretion profile (||) 

A{r) =A(r Q ) e'^f- 

can be rescaled so that at the end of the simulation 
a(rQ,tc) = 50 Mq pc~ 2 . This zero-point does not in- 
fluence the profile, nor the chemical evolution, so it can 
always be rescaled a posteriori. 
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For our example, we take as the reference static case a 
model adopting a Schmidt SF law with k = 1.5 (Kennicutt 
1998) and a uniform infall time-scale of 3 Gyr (model S15a 
of PC99, see also Table @). 

4-.1. Models with inflow from the outer disc 

In model S15RFa a uniform radial inflow pattern with 
v = —1 km sec -1 and inflow from the outer disc is im- 
posed upon the static model S15a, with no further change 
in the model parameters (Table [l]). Fig. shows the ef- 
fects of inflow on the total surface density and gas density 
distribution, and on the oxygen gradient (details on the 
observational data in the plots can be found in PC99). By 
comparing the solid line and the dashed line, we notice 
the following main effects. 

— Since matter flows inward and accumulates toward the 
Galactic Centre, the density profile gets steeper in the 
inner parts, while in the outer parts it remains rather 
flat because gas is continuously poured in by the flat 



outer gaseous disc. The gas density distribution shows 
a similar behaviour. 
— The overall metallicity gets much lower because of the 
dilution by primordial gas inflowing from the outer 
disc. The gradient becomes steeper, especially in the 
outermost shells, since the discontinuity in metallicity 
between the star-forming disc and the outer gaseous 
layer is smeared inward by radial inflows. 

To compensate for the steepening of the density distribu- 
tion induced by radial inflows, we must adopt a shallower 
initial accretion profile; our simulations show that a scale 
length rd ~ 5 kpc for the infall profile reduces in the end 
to a density profile matching the desired scale length of 
~4 kpc at the Solar Neighbourhood. At the same time, the 
chemical enrichment must get more efficient for the overall 
metallicity to increase to the observed levels; the predicted 
metallicity is improved by adopting an IMF more weighted 
towards massive stars, i.e. by increasing the "IMF scaling 
fraction" £ (see PCB98 and PC99 for a description of our 
model parameters). Together with the SF efficiency v, C 
is re-calibrated to match the observed gas surface density 
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and metallicity at the Solar Neighbourhood (as for the 
models of PC99). In this way we calibrate model S15RFb 
with respect to the Solar Neighbourhood; see Table g for 
details on the adopted parameters. Comparing now this 
re-calibrated model with radial flows (dotted line) to the 
original model S15a, the metallicity gradient is increased 
with respect to the static case, but still a bit flat with re- 
spect to observations. The gas density distribution peaks 
in the inner regions, as expected, and remains quite high 
(much higher than observed) in the outer regions due to 
substantial replenishment from the outer disc. 

4-2. Models with no inflow from the outer disc 

Let's now consider the case when radial flows are limited 
within the star forming disc and there is no inflow from 
the external gaseous disc (v N+ i = 0). If radial flows are 
mainly driven by shocks in spiral arms, for instance, they 
might indeed occur only within the stellar disc, while the 
outer gas layer remains unperturbed. 

Fig. H (left panels) shows how the final density profile 
becomes much steeper than the reference accretion profile 
(model S15RFc, dashed line), as expected since matter 
is efficiently drifting inward with no replenishment from 
the outer disc. For the final density profile to match the 
observed one, we must assume a much shallower initial 
accretion profile (r^ ^6 — 7 kpc) . Then the resulting local 
density profile is close to the observed one, while in the 
outer parts the profile remains steeper (model S15RFd, 
dotted line). 

The gas density profile shows a similar behaviour, 
strongly peaked toward the centre while dropping quickly 
(much more quickly than observed) outside the Solar ring 
(upper right panel). 

The overall metallicity is reduced with respect to the 
original model S15a (lower right panel, model S15RFc) 
and again we need to increase the IMF scaling fraction £ 
to raise the chemical enrichment to the observed values 
(model S15RFd). The resulting gradient is roughly com- 
parable to the observed one. 

4-3. Concluding remarks 

We can summarize the general effects of radial flows as 
follows. 

— Since matter flows inward and accumulates toward the 
Galactic Centre, the density profile in the inner parts 
gets steeper. A shallower intrinsic accretion profile is 
to be adopted, in order to recover the observed local 
scale-length at the end of the simulation. 

— In the outer parts, the profile declines more or less 
sharply, depending on whether the inflow occurs only 
within the star-forming disc or there is also substantial 
inflow from the outer, purely gaseous disc. 



— The gas distribution shows a similar behaviour: it re- 
mains quite flat in the outer regions if gas is poured 
in by the flat outer gas disc, while it tends to drop 
sharply (more than observed) otherwise. 

— The inner gas profile is very steep in the case of a 
Schmidt SF law (models S15RF) — while with other 
SF laws with radially decaying efficiency (see PC99 
and §||) the effect would be somewhat compensated by 
a larger gas consumption by SF in the inner regions. 

— The overall metallicity gets much lower because of the 
dilution by gas inflowing from metal poor outer shells 
(and possibly the primordial outer disc). A higher frac- 
tion of stars contributing to the chemical enrichment is 
needed in the model to match the observed metallicity 

— The metallicity gradient tends to steepen, in agree- 
ment with results by other authors (see references in 

i). 

5. Some "successful" models 

In §|] we have illustrated the qualitative effects of ra- 
dial inflows using, for the sake of example, models with a 
Schmidt SF law. We will now consider models with various 
SF laws (see PC99) and tune the inflow velocity pattern, 
in each case, so as to match the observational data on 
the radial profile of the Galactic Disc. The various "suc- 
cessful" models presented here should not be taken as de- 
tailed, unique recipes to reproduce the Disc. Rather, they 
are meant as examples of how inclusion of radial flows in 
the chemical model can improve the match with the data. 

5.1. Models with Schmidt law 

Let's first consider again the Schmidt SF law. Inspection 
of models S15RFb and S15RFd (with and without radial 
inflows from the outer disc) suggests to proceed as follows. 



S15RFe 




-0.2 - 
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R 

Fig. 4. Inflow velocity pattern for the "successful" models 
S15RFe and OlORFe. 
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Fig. 5. A "successful" S15RFe model compared to the case Fig. 6. A "successful" O10RF model compared to the case 
with no flows. Data and symbols as in Fig. 0. with no flows. 



— Some inflow from the outer gaseous layer is needed to 
reproduce the shallow decline of the gas distribution 
observed out of the Solar ring, though the inflow speed 
from the outer disc should be slower than -1 km sec -1 
otherwise the predicted gas density is too high in the 
outer parts (as in model S15RFb). 

— Models with radial inflows can predict metallicity 
gradients close to the observed ones even with a 
Schmidt SF law, provided drift velocities within the 
star-forming edge are relatively high (of the order of 
-1 km sec" 1 ). Inflow patterns decelerating inward are 
particularly efficient in building the metallicity gradi- 
ent (Gotz & Koppen 1992). 

An example of a drift velocity profile shaped according to 
these considerations is shown in Fig. [|; the corresponding 
model gives indeed a reasonable fit to the data (model 
S15RFe, Fig. ||). The gas profile keeps increasing inward, 
yet the model does not reproduce the peak of the molec- 
ular ring, for which we need to include the peculiar radial 
flows induced by the Bar (see the discussion in PC99 and 
§||). Model S15RFe shows how radial inflows can in fact 
allow for metallicity gradients comparable with the ob- 



served ones, even with a Schmidt SF law which would be 
excluded on the sole base of static models (PC99). 

5.2. Models with spiral-triggered SF laws 

Let's consider now models adopting an Oort-type SF law, 
with a SF efficiency inversely proportional to the galacto- 
centric radius and a Schmidt-like exponent k = 1.0 (Ken- 
nicutt 1998, PC99). The static and the "successful" model 
with this SF law are model OlOa from PC99 and model 
OlORFe, respectively (Table 0). 

As in the static case, (see PC99), these models behave 
somewhat similarly to models adopting the Schmidt SF 
law. Fig. [j] shows that a good fit is obtained with model 
OlORFe, where we applied the velocity pattern of Fig. |^ 
with inflows becoming slower inwar, very close to that 
used for model S15RFe. Notice that in this case the higher 
SF efficiency in the inner region and the slowdown of in- 
flows conspire to accumulate gas around r = 3 kpc while 
consuming it at inner radii, generating a peak in the pre- 
dicted gas distribution which closely reminds the observed 
molecular ring. 
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Fig. 7. Inflow velocity pattern for the "successful" model 
DRRF 



5.3. Models with gravitational self-regulating SF laws 

We now investigate models adopting a self-regulating SF 
process driven by gravitational settling and feed-back 
from massive stars, implying a SF efficiency exponentially 
decaying outward in radius, such as the law by Dopita & 
Ryder (1994; see PC99). The reference static model in this 
case is model DRa of PC99, while the corresponding "suc- 
cessful" model with radial flows is DRRFe (see Table [l]for 
the relevant parameters). 

To reproduce the observed metallicity gradient with 
this SF law, negligible inflow is required in the outer re- 
gions where static models already predict a gradient with 
the right slope (see model DRa), while a moderate inflow 
velocity (^0.3 km sec -1 ) in the inner parts is needed to 
maintain the observed slope also where the predicted gra- 
dient would otherwise flatten (see PC99). Such a model is, 
for example, model DRRFe, obtained with the inflow ve- 
locity pattern shown in Fig. (?]. In Fig. || the "successful" 
model DRRFe is compared to the original model DRa of 
PC99 with no radial flows. 



5.4- Concluding remarks 

The aim of the previous "successful" examples is to remark 
how, for each SF law considered, it is possible to tune 
the inflow velocity pattern so as to get a good overall 
agreement with the observational data. As an indication, 
suitable inflow profiles have the following characteristics. 

— With a Schmidt SF law, relatively fast (but still plau- 
sible) flows are needed within the star-forming disc, 
with velocities of the order of 1 km sec -1 . The required 
metallicity gradient is easily obtained, especially if the 
inflow velocity is slowly decreasing inward. 
Inflow from the outer gaseous disc must be less promi- 
nent (though present), otherwise the predicted gas 
density in the outer regions is too high. 
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Fig. 8. A "successful" DRRF model compared to the case 
with no flows. 



— In models adopting an Oort-type SF law with k = 1.0, 
a good fit to the data is obtained with an inflow 
profile of the same kind: inflow at ~ 0.5 km sec -1 
from the outer gaseous disc, drift velocities raising at 
~ 1 km sec -1 at the stellar disc edge and declining 
inward. Curiously enough, in this case the higher SF 
efficiency in the inner regions combines with the slow- 
down of inflows to predict a peak in the gas distribu- 
tion around r = 3 kpc, quite reminiscent of the ob- 
served molecular ring. 

— Models adopting the SF law by Dopita & Ryder give 
a good description of the outer regions of the Disc al- 
ready in the absence of flows (PC99). Mild radial in- 
flows can be assumed in the inner regions, where the 
radial gradient would otherwise flatten, to obtain the 
right slope throughout the disc. The required drift ve- 
locities are around 0.3 km sec -1 , which can be reason- 
ably provided, for instance, by the shocks occurring 
within the spiral arms (see §^). 

The various "successful models" presented here are not 
meant to be definitive recipes to reproduce the radial pro- 
file of the Galactic Disc. In fact, the adopted inflow veloc- 
ity profiles are quite arbitrary and ad hoc. These models 
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are rather meant to show how radial flows can be a viable 
mechanism to interpret the properties of the Disc. In PC99 
we showed how none of the various SF laws investigated is 
able, by itself, to reproduce the observed metallicity gradi- 
ent throughout the whole extent of the Disc, unless some 
additional "dynamical" assumption is included. In PC99 
we considered the classical case of an inside-out disc for- 
mation, namely of an infall time-scale increasing outward. 
Here, we just show that radial inflows can provide another 
viable "dynamical" assumption to be combined with any 
SF law to reproduce the observed gradient. In particular, if 
we adopt a Schmidt SF law with k = 1.5 or an Oort-type 
SF law with k = 1.0, as recent empirical evidence seems 
to support (Kennicutt 1998), the required variation of the 
accretion time-scale is too extreme in the pure inside-out 
assumption (PC99), and radial inflows are then necessary 
to explain the metallicity gradient. Of course, the two ef- 
fects (inside-out formation and radial flows) can also play 
at the same time. We do not address here their combined 
outcome, because in our models this would merely trans- 
late into increasing the number of parameters which one 
can tune to fit the observational data. No further insight in 
the problem would be gained. With this kind of models we 
can just learn how the various players (different SF laws, 
radially varying accretion time-scales, radial gas inflows) 
enter the game of reconstructing the general picture, and 
analyse their effects one by one. 

It is also worth stressing that even slow radial gas 
flows, with velocities well plausible in terms of the trigger- 
ing physical mechanisms and within the observational lim- 
its (§^|), have non-negligible influence on chemical models, 
especially on the gas density distribution. It is therefore 
misleading to seek for a one-to-one relation between gas 
content and metallicity, or between gas profile and metal- 
licity gradient, like the one predicted by simple models 
(e.g. Tinsley 1980). When studying the chemical features 
of galaxies, it is very dangerous to assume that such a re- 
lation must hold, since even mild flows can easily alter the 
overall distribution. 

The models with smooth radial inflows presented in 
this section (with the possible exception of model OlORFe) 
are still unable to reproduce the gas density peak corre- 
sponding to the molecular ring around 4 kpc, since that 
needs to take into account the peculiar dynamical influ- 
ence of the Galactic Bar on gas flows. This will be the 
issue of the next section. 

6. The role of the Galactic Bar 

There is by now substantial evidence that the Milky Way 
hosts a small Bar in its inner 3 kpc or so. The idea 
was originally suggested to explain the kinematics of the 
atomic and molecular gas near the Galactic Centre (de 
Vaucouleurs 1964, Peters 1975, Liszt & Burton 1980, Mul- 
der & Liem 1986). In recent years further evidence for 
a Galactic Bar has piled up from several tracers: gas 



dynamics (Binney et al. 1991; Weiner & Sellwood 1996, 
1999; Yuan 1993; Wada et al. 1994; Englmaier & Ger- 
hard 1999; Fux 1999), IR photometry and star counts 
(Blitz & Spergel 1991, Weinberg 1992, Nikolaev & Wein- 
berg 1997, Dwek et al. 1995, Binney et al. 1997, Unavane 
& Gilmore 1998), stellar kinematics (Zhao et al. 1994; 
Weinberg 1994; Fux 1997; Sevenster 1997, 1999; Seven- 
ster et al. 1999; Raboud et al. 1998), OGLE data (Stanek 
et al. 1994, 1997; Paczynski et al. 1994; Evans 1994; Zhao 
et al. 1995, 1996; Zhao & De Zeeuw 1998, Ng et al. 1996). 
For a review on the Galactic Bar see e.g. Gerhard (1996, 
1999). The determination of the characteristic parameters 
(size, axis ratio, rotation speed, orientation and so forth) 
is even more difficult for the Bar of our own Milky Way 
than for external galaxies. Broadly speaking, the various 
studies mentioned above indicate a Bar with a major axis 
of 2-4 kpc viewed at an angle of 15-45° in the first lon- 
gitude quadrant, an axis ratio around 3:1 and a pattern 
speed ftp ~ 60 km sec -1 kpc -1 . 

PC99 underlined that the dynamical influence of the 
Galactic Bar is likely to account for the peak at 4 kpc dis- 
played by the gas profile in the Disc, which static chemical 
models are unable to reproduce if other constraints, like 
the observed metallicity gradient, are to be matched as 
well. In fact, gravitational torques in a barred, or non- 
axisymmetric, potential are thought to induce gas ac- 
cumulation and formation of rings at the corresponding 
Lindblad resonances (e.g. Combes & Gerin 1985; Schwarz 
1981, 1984). In brief, bar-induced flows sweep gas away 
from the co-rotation (CR) radius, where the bar roughly 
ends, toward the inner and outer Lindblad resonances 
(ILR, OLR). In fact, we developed our new chemical model 
with radial flows also with the aim to mimic such effects of 
the Bar upon the gas distribution, by simulating suitable 
flow velocity profiles. 

As mentioned above, though the existence and gross 
features of the Galactic Bar are by now established, there 
is no general agreement on details like its size and pat- 
tern speed, and on the corresponding radii for its CR and 
ILR, OLR. In this paper, with the aim to reproduce the 
molecular ring at 4-6 kpc, we will consider two Bar models 
covering the range of scenarios suggested in literature: 

case A) Bar's CR around 3.5 kpc, so that the dip in the 
gas distribution between 1.5 and 3.5 kpc is interpreted 
as a fast inward drift of gas from CR to the ILR and 
the nuclear ring; 

case B) Bar's CR around 2.5 kpc and OLR around 
4.5 kpc, so that the molecular ring is interpreted as 
accumulation of gas from CR toward the OLR. 

In any case, we will assume here that the Bar influences 
only the inner 5-6 kpc of the Galactic Disc, where its OLR 
is supposed to lie at the outermost according to current 
understanding, while leaving regions outside the OLR un- 
affected (see also Gerhard 1999). Actually, the possible 
influence of the Bar over a larger Disc region than its for- 
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T b»r= 5 
T b,r= 1 

no Bar 



models S15RF 
Bar A 



Fig. 9. Inflow velocity pattern adopted from the onset of 
the Bar at the time to ~ 2b or upon the base "no Bar" 
model (S15RFe), to mimic case A for the Galactic Bar 
with models S15RF. The dotted line applies to the base 
model S15RFe, and also to the "barred" models before 

to — Tb ar . 

mal extent is still an open problem, as we will comment 
upon in the final conclusions (§0). 

In the framework of our chemical model, the inclusion 
of the Bar translates into imposing a suitable velocity pro- 
file for the radial flows in the inner regions of the disc. 
Namely, we run the "successful models" with radial flows 
presented in §[|, but at a suitable age the Bar is assumed to 
form and the radial velocity profile is altered correspond- 
ingly, by modifying the coefficients (3k, 7& describing 
the radial flow pattern (§||). In case A, we will impose fast 
radial inflow velocities within CR at 3-4 kpc to mimic the 



rapid drift of the gas toward the ILR (£6.1). In case B, we 



w ill i mpose outflows from CR to the OLR around 4.5 kpc 

As to the age of the Bar, an upper limit is set by 
the typical age of its stellar population, 8-9 Gyr (Ng 
et al. 1996); an age of 5 to 8 Gyr has been suggested by 
Sevenster (1997, 1999). The results from our simulations 
turned out to be quite insensitive to a change in the Bar's 
age from 9 to 5 Gyr; therefore, we will present simulations 
with a 5 Gyr old Bar (Tb ar = 5) as representative of the 
generic case of age > 5 Gyr. We will also consider, for the 
sake of completeness, the case of a much younger Bar of 
1 Gyr of age (T bar = 1). 

6.1. Modelling the effects of the Bar: case A 

To mimic case A, where the Bar induces fast inflows from 
CR to its ILR, the inflow velocity is typically increased for 
r < 3.5 — 4 kpc (the Bar's CR radius) with respect to the 
drift pattern adopted before the onset of the Bar. Mod- 
els corresponding to case A give good results when com- 
bined with a Schmidt SF law with radial inflows. Starting 
from the corresponding "successful" model S15RFe of §|j 
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Fig. 10. Models S15RF mimicking case A for the Galactic 
Bar, adopting the inflow patterns of Fig. ^|. The "no Bar" 
model S15RFe is also shown for comparison as a dotted 
line. 

and changing the velocity law as in Fig. [9] at the onset of 
the Bar, we obtained the models shown in Fig. |l(| com- 
pared to the original model S15RFe with no Bar's effects 
included. Obviously, if the Bar is younger (lf, or = 1 Gyr) 
faster induced inflows need to be assumed in the inner re- 
gions so that the observed sharp dip in the gas profile at 
r < 3.5 kpc is obtained in a shorter time (Fig. ||). Any- 
ways, no extreme speeds need to be induced by the Bar 
(\v\ < 5 km sec) at these radii yet, to resemble the ob- 
served gas profile, so the models remain plausible. 

This type of solution actually corresponds to the one 
originally suggested by Lacey & Fall (1985), who in fact 
assumed a Schmidt SF law, radial inflows of the order of 
-1 km sec -1 to reproduce the metallicity gradient, and for 
r < 4 kpc a raise in the inflow speed up to -10 km sec" 1 to 
reproduce the gas profile. Without such a spike in the in- 
flow velocity, the gas distribution keeps rising inward, with 
no depression (cfr. model S15RFe). This picture also corre- 
sponds to the situation for viscous models as suggested by 
Thon & Meusinger (1998): the detailed gas profile of the 
Disc can be reproduced only by artificially increasing the 
viscosity in the inner regions, so as to mimic the influence 
of the Galactic Bar. 

Case A can also well combine with model CTORFe, 
namely with an Oort-type SF law with k = 1 plus the 
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Fig. 11. Inflow velocity pattern adopted from to — Tbar 
upon the base ("no Bar") model OlORFe, to mimic case 
A for the Galactic Bar with models O10RF. 
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inside 4 kpc one can fit the observed gas depression. As an 
example, we show models O10RF with inflow patterns as 
in Fig. [ll]; the corresponding results are plotted in Fig. [l^ 
together with the base model OlORFe. 

Notice how in all these "case A" models the metallicity 
gradient is only negligibly affected by the switch of the in- 
ner inflow pattern at the time of onset of the Bar, t G — Tbar 
(to — 15 Gyr is the present Galactic age). The effects are 
at most limited to r < 3 kpc, where they are hard to check 
from observations, since in that region the gas distribu- 
tion is depressed and the metallicity tracers are missing 
as well, since they are young objects strongly correlated 
to present-day SF activity and therefore to the presence 
of gas. We remark that abundance data at r = in the 
plots for the abundance gradient are to be disregarded as 
a constraint for the model, since they refer to the Galactic 
Centre population, not to the Disc (see PC99). 

6.2. Modelling the effects of the Bar: case B 

According to what we labelled above as "case B", the 
Galactic Bar ends in correspondence to its CR around 
2.5 kpc, and it has an OLR around 4.5 kpc. Therefore, the 
gas is expected to drift from CR outward and accumulate 
at the OLR, while gas inflowing from outer regions slows 
down and accumulates as well at the level of the resonance 
at 4.5 kpc. 

Models for case B will adopt, starting from t G — T bar , a 
radial flow pattern with positive velocities (outflow) from 
r = 2.5 kpc to r ~ 4.5 kpc, and negative velocities (inflow) 
from outside dropping to zero at r ~ 4.5 kpc. Case B can 
be combined with model DRRFe, namely with the SF law 
by Dopita & Ryder (1994) and the corresponding "success- 



Fig. 12. Models O10RF mimicking case A for the Galactic 
Bar, adopting the inflow patterns of Fig. [II]. 




corresponding "successful" inflow pattern. In fact, also in 
model OlORFe (which is actually the only model able to 
predict a peak in the gas profile reminding of the molec- 
ular ring even without including the Bar, see §[s]) the gas 
profile rises inward, and by increasing the inflow velocities 



Fig. 13. Radial flow velocity pattern adopted from t G — 
Tbar upon the base "no Bar" inflow pattern, to mimic case 
B for the Galactic Bar with models DRRF. The dot-dashed 
horizontal line at v = marks the transition from inflows 
(negative v) to outflows (positive v). 



14 



Portinari and Chiosi: On radial gas flows, the Galactic Bar and chemical evolution in the Galactic Disc 



models DRRF 
Bar B 




Fig. 14. Models DRRF mimicking case B for the Galac- 
tic Bar, adopting the flow patterns of Fig. |l3|. Disregard 
abundance data at r = (see text). 



ful" overall inflow pattern (§g). The relevant models with 
their detailed velocity patterns, for the usual two values 
of Tbar, are shown in Figs. [l3| and O. Notice once again 
that rather low drift velocities (|u| < 0.5 km sec) suffice to 
reproduce the gas peak, which reinforces the plausibility 
of the models. 

In case B models, the metallicity gradient is clearly af- 
fected in the inner regions by the switch in the gas flow 
pattern, at least if this lasted for some time (Fig. [l4|, case 
Tbar = 5). This feature, though, is again hard to check 
from observations since there are no tracers of Disc metal- 
licity for r < 4 kpc (see 



6.1) 



Case B can also reproduce the observed gas profile 
combined with an Oort-type SF law with k = 1 (model 
OlORFe), especially provided the Bar formed recently, as 
displayed in Fig. |l^ with Tb ar = 1. If the Bar-induced 
flows activate much earlier, the predicted peak is very 
sharp and narrow (Fig. (l6L case Tbar- — 5 Gyr), due to 
the milder sensitivity of SF to the gas density in this case 
(Schmidt-like exponent k — 1): with respect to other SF 
laws, this SF is relatively less efficient where the gas den- 
sity is high, namely where the gas accumulates around 
4.5 kpc, while it is relatively more efficient where the den- 
sity drops, below 4 kpc. Such a sharp peak seems in con- 




Fig. 15. Radial flow velocity pattern adopted from to — 
Tbar upon the base "no Bar" inflow pattern, to mimic case 
B for the Galactic Bar with models O10RF. 
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Fig. 16. Models O10RF mimicking case B for the Galactic 
Bar, adopting the flow patterns of Fig. 15. 



trast with the observed, quite broad distribution (the ob- 
servational uncertainty on the gas density profile in the in- 
ner Galactic region is less than a factor of 2, Dame 1993). 
Models O10RF with an "old" Bar are therefore less ap- 
pealing in case B. 
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6.3. Concluding remarks 

We referred to the current understanding of the structure 
and features of the Galactic Bar to simulate its dynamical 
influence on the surrounding gas flows and distribution. 
Chemical models accounting for the effect of the Bar make 
it finally possible to reproduce the gas profile properly, 
which could not be accomplished by static models (PC99). 
Broadly speaking, two main scenarios are presented. 

1. The Bar stretches to 3.5 kpc and the gas peak is due 
to a rapid depletion of gas drifting from the Bar's CR 
to the Galactic Centre (case A). To reproduce the gas 
profile, we need a gas distribution which keeps increas- 
ing from the outer disc inward, down to ~ 4 kpc where 
the Bar's influence produces the sharp depletion. This 
can be achieved only by means of efficient radial in- 
flows over the whole disc, and in this case radial inflows 
are also the major source for the observed metallicity 
gradients (models S15RFe and OlORFe). 

2. The Bar ends around 2.5 kpc, where its CR is set, and 
the gas peak is due to the accumulation of gas from CR 
outward to its OLR at ~ 4.5 kpc (case B). A rather 
limited contribution of radial inflows from the outer 
regions suffices to obtain the peak. In this case, the 
metallicity gradients in the outer regions of the disc 
may be mainly due to the SF process itself, and to 
the intrinsic variation of the SF efficiency with radius 
(model DRRFe). 

Within these simplified models it is unfeasible to dis- 
cuss any further on the scenarios for Bar structure and 
age, and related gas flows. Only detailed dynamical simu- 
lations for Bar formation, evolution and potential can tell 
how the molecular ring consequently formed, which is be- 
yond the goals of this paper (see also §0) . Here we were just 
interested in showing how even simple qualitative models 
for Bar-induced gas flows in the inner disc solve in fact the 
puzzle encountered in PC99. Namely, they can reproduce 
at the same time both the metallicity gradient and the gas 
distribution, in particular the peak corresponding to the 
molecular ring around 4 kpc. This can be accomplished 
already with quite slow, and largely plausible, flow veloc- 
ities. The present modelling therefore provides a simple 
tool for qualitative understanding of possible behaviours. 

Regardless of details, however, one condition is neces- 
sary for the scheme to work: there must be enough gas 
in the inner regions of the Disc, which the Bar can then 
"shape" to resemble the detailed observed distribution. 
This favours chemical models with radial inflows where 
the metallicity gradient can coexist with high gas frac- 
tions in the inner shells. Our simulations showed that no 
Bar-induced gas flow superimposed on otherwise static 
models (as those by PC99) can produce the observed gas 
peak: whatever the assumed age of the Bar or gas velocity 
profile, there is not enough gas left in the inner Galac- 
tic regions if we are to reproduce the metallicity gradient 



as well. Gas must be continuously replenished by inflows 
from outer regions; that's why we presented here "barred" 
models based only on the "successful" models with radial 
inflows from §|| and never on the static models of PC99. 

7. Summary and conclusions 

From the results of static chemical models, PC99 under- 
lined the need to introduce radial flows to explain some 
features of the Galactic Disc. In fact, static models are un- 
able to reproduce, at the same time, both the metallicity 
gradient and the radial gas profile; in particular, the peak 
corresponding to the molecular ring at 4-6 kpc is likely to 
be a consequence of gas drifts induced by the dynamical 
influence of the Galactic Bar. Therefore, in the present pa- 
per we introduced a new chemical model including radial 
gas flows, developed as a multi-dimensional generalization 
of the original static model (§|3|). Our model is conceived so 
as to adapt to any imposed radial velocity profile, describ- 
ing both inflows and outflows in any part of the disc. The 
model is carefully tested against instability problems and 
spatial resolution, by comparing it to suitable exact ana- 
lytical cases (Appendix B). In this paper we applied the 
model to the Galactic Disc; more in general, such mod- 
els allowing for gas drifts are meant to be used as fast 
and handy interfaces between detailed dynamical galaxy 
models (predicting the velocity profiles) and parametric 
chemical and spectro-photometric models. 

An overview of the behaviour of chemical models with 
radial inflows of gas shows that these provide an alterna- 
tive "dynamical" assumption to the inside-out disc forma- 
tion scenario to explain the metallicity gradient (§|J). With 
radial gas flows, the model can reproduce the metallicity 
gradient even in the case of a Schmidt or an Oort-type 
SF law, which were excluded in the case of static models 
(see PC99). In addition, it appears that even low radial 
flow velocities, well within observational limits and the- 
oretical expectations (see §j|), have non-negligible effects 
upon model predictions on the metallicity gradient and 
moreover on the gas distribution. In particular, if radial 
gas inflows are allowed for, a metallicity gradient can co- 
exist with a high gas fraction in the inner regions, at odds 
with simple static models. This is indispensable to repro- 
duce the observed gas distribution in the inner Galaxy 
(see point 2 below). The remarkable effects of even slow 
radial flows upon observable quantities, mainly upon the 
gas distribution, should be kept in mind as a caveat when 
comparing real galaxies to simple analytical models which 
predict a one-to-one relationship between metallicity and 
gas fraction (e.g. Tinsley 1980). Our models show that 
small dynamical effects, like slow gas flows, can easily 
make real systems depart from the behaviour of simple 
models. 

With our model it is possible to mimic the dynamical 
influence of the Galactic Bar and reproduce the peak in 
the gas distribution around 4 kpc (§§). Two scenarios, 
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related to two different models for the structure of the 
Bar, are qualitatively suggested. 

A. With a Schmidt or an Oort-type SF law, slow ra- 
dial inflows in the disc pile up gas inward down to 
r = 3.5 — 4 kpc. Here, the Bar CR radius is found and 
the gas is quickly swept inward from CR toward an 
ILR, which causes the drop in the gas profile at 3.5 kpc. 

B. With a SF law like that by Dopita & Ryder (1994), 
smaller inflow rates suffice to reproduce the metallicity 
gradient, leading to a lower concentration of gas in 
the inner regions than in the previous case. The peak 
corresponding to the molecular ring can be reproduced 
with a Bar CR around 2.5 kpc and its OLR around 
4.5 kpc, so that all the gas external to r — 2.5 kpc 
tends to pile up around the OLR. 

Though these models are just qualitative and cannot de- 
scribe the detailed dynamical process of Bar formation nor 
the evolution of the related gas flows to form the molecular 
ring, they provide two interesting indications. 

1. Only when introducing the effects of the Bar, the 
model is able to reproduce the radial gas profile prop- 
erly. The only possible exception resides in a particular 
combination of an Oort-type SF law with a radial in- 
flow pattern whose velocity decreases inward (model 
OlORFe; this combination may lead to a peak of the 
gas distribution in the inner Galactic regions, closed to 
the observed molecular ring. But this particular, for- 
tunate case does not diminish the general conclusions 
about the role of the Galactic Bar. 

2. In any case (A or B above), overall radial inflows in the 
disc are indispensable to replenish the inner regions 
with enough gas that the observed molecular ring can 
form under the influence of the Bar. This seems to 
favour disc models with radial inflows, unless one as- 
sumes that the gas in the ring has some different origin 
(gas swept from the Bulge, or accreted later). 

To investigate these issues any further, detailed gas- 
dynamical models are obviously required. Unfortunately, 
most studies on Bar-induced gas dynamics (see references 
in §||) concentrate on the observed features of the very in- 
ner regions, such as the nuclear ring, the 3 kpc expanding 
arm, and so forth. Little discussion can be found about 
the effects of the Bar on more external regions, and on 
the formation of the molecular ring in particular: whether 
it is due to gas depletion inside CR as in our case A, or 
due to gas accumulation at some resonance (e.g. Binney 
et al. 1991, Fux 1999) as in our case B, or whether it just 
consists of two or more tightly wound spiral arms (e.g. 
Englmaier & Gerhard 1999). Further gas-dynamical stud- 
ies suggesting detailed scenarios, time-scales, and velocity 
profiles for the formation of the molecular ring would be 
welcome, for the sake of including the effects of the Bar in 
chemical evolution models more consistently. 



Further investigation of gas-dynamical models on the 
influence of the Bar on even larger scales (namely, outside 
its OLR) should be pursued as well, since this is a clue 
issue related to a claimed discrepancy between the char- 
acteristics of the Galactic Bar and the observed metallicity 
gradient. It is well known that barred galaxies display sys- 
tematically shallower gradients than ordinary spirals (e.g. 
Alloin et al. 1981, Vila-Costas & Edmunds 1992, Martin 
& Roy 1994). This is likely a consequence of the radial mix- 
ing induced by bars; in fact, Martin & Roy (1994) found 
a correlation for external galaxies between the strength 
of a bar and the metallicity gradient. Taking this empir- 
ical relation at face value, the Galactic Bar with an ax- 
ial ratio of ~ 0.5 should induce a metallicity gradient of 
-0.03 dex/kpc, much shallower than the observed one of 
-0.07 dex/kpc, which is typical of a normal Sbc galaxy. 
To overcome such a puzzle, it has been suggested that the 
Galactic Bar must be very young (<1 Gyr), so that there 
was not enough time yet to flatten the gradient (Gum- 
mersbach et al. 1998); but this is in conflict with other 
estimates of the Bar's age (e.g. Sevenster 1997, 1999). Al- 
ternatively, we suggest that the discrepancy might be only 
apparent, since the Galactic Bar is quite small, and the 
Milky Way cannot be properly considered a barred spiral. 
It might be unlikely that the Bar can influence the metal- 
licity gradient all over the Disc, as in really barred galax- 
ies: Bar-induced radial drifts and corresponding chemical 
mixing are expected to occur from CR toward the ILR (in- 
flows) and to the OLR (outflows; e.g. Schwarz 1981, 1984; 
Friedli et al. 1994). Present understanding of the Galactic 
Bar sets its OLR between 4.5 and 6 kpc (§|6]and references 
therein), so in our models we presumed that the Bar in- 
duces negligible mixing beyond these radii, regardless of 
its age (see also Gerhard 1999). If the situation is as in the 
models we presented here, in fact, the metallicity gradi- 
ent in the outer regions is unperturbed and just related to 
intrinsic Disc properties and/or large-scale viscous flows. 

However, gas-dynamical simulations dedicated to the 
effects of the Galactic Bar over the whole Disc would 
be necessary, so as to investigate the relation between 
the Bar, radial mixing and the metallicity gradient more 
consistently. More in general, including the effects of 
bar-induced radial flows in the picture of the chemical 
evolution of spiral galaxies might turn out to be of wide 
interest, since it is likely that all spirals develop at some 
point, or have developed in the past, some bar-like struc- 
ture (Binney 1995). Infrared observations indeed reveal 
that a large fraction of spirals host a barred structure (e.g. 
Eskridge et al. 1999), and recent numerical simulations 
suggest that even weak bars or oval distortions may be 
able to induce radial drifts to form multiple gaseous rings 
at the corresponding Lindblad resonances (Jungwiert & 
Palous 1996). Bars could even drive secular evolution of 
spiral discs from late to early type (e.g. Dutil & Roy 1999). 
Bar-driven radial gas flows might therefore play a funda- 
mental role in the chemical evolution of spiral discs. 



17 



Acknowledgements. We thank Joachim Koppen for his advice 
on numerical modelling of radial gas flows, Yuen K. Ng and 
Antonella Vallenari for useful discussions about Galactic struc- 
ture, and our referee, Mike Edmunds, whose suggestions much 
improved the presentation of our paper. 

L.P. acknowledges kind hospitality from the Nordita Institute 
in Copenhagen, from the Observatory of Helsinki and from 
Sissa/Isas in Trieste. This study has been financed by the 
Italian MURST through a PhD grant and the contract "For- 
mazione ed evoluzione delle galassie", n. 9802192401. 

Appendix A: The explicit expression of Eq. ( |Tg| ) 

Here we sort out the explicit expression ( po| ) of the solution (^) 
for Gi(rk,t), by calculating the matrix e tA . Let us first intro- 
duce the simplified notation: 



where we have indicated with: 



G h {t) = Gi(r h ,t), 



W h = Wi(r k ) 



(Al) 



For the sake of example, we present here the detailed solution 
in the case of three shells (TV = 3). The relevant system of 
equations drij ) becomes: 

f t Gi(t) = tfiGi(t) + 7 iG 3 (t) + Wi 

£ G 2 (t) = a 2 Gi(t) + 2 G 2 (t) + 72G3W + W 2 (A2) 

^ Ga(t) = as G 2 it) + ^G 3 it) + W 3 +— l 

where all the coefficients are considered as constants (see § [j.2j ) , 
though we hav e om itted the bar over Wt and for simplicity. 
A system like (A2) holds for each chemical speci es i, but the 
characteristic matrix A is the same for any i (see £3.2); in this 
case: 

/ 0i 71 \ 
A = I Ul -&2 72 

\ « 3 t) 3 / 

Let's first notice, from the definition (^|), that 7^ and ctk+i can 
never be both positive: at least one of them must be zero since 
they are "activated" in the opposite cases of inflow or outflow 

+ 5' 



at r k+ i, respectively; if there is no flow at all through r k 
they both reduce to zero. Therefore, in our calculations we are 
always entitled to use the condition: 



Jk ctk+i =0 Vfc 

The eigenvalues of the matrix A are 

Al,2,3 = 01,2,3 

and the associated eigenvectors are: 
/ 1 

Ul = 



(A3) 



/ 71 

/ *2-#l 



c*2 



^3 a 2 



u 2 = 



From @ we get: 

02/21 (*) 



1 



U 3 = 



/ 21 72 , 

/ 03—0! 3 -0 2 \ 

/ 



72 



# 3 -#2 
1 



71/21W s^[hi(t) - hi(t)Y 

e fl2t 72/32(t) 



\^^[/32(t)-/ 3 i(i)] oa/saW 



fklit) = flkit) 

Defining: 



0fc-0! 



'to 



»«,(*!-*) 



I?* 



the solution (|l^) in the case N — 3 becomes: 

r Gi(ti) 



e" lAt Gi(t ) + 



+ 71 



„tfo At *9i At 



■G 2 (t ) + 



+ 7i72 ^%^ G 3 (to) + 



G? a (*i 



+ ^-V, W * [ g 3 (At)- 91 



Q 2 <J-1 1*0 ) + 

+ e^ At G 2 (* ) + 



+ 



g3 (At)- gl (At) _ 32 (At)- gi (Ai)1 



9 3 At_ e >),At 



+ 72 ' 

+ a 2 



+ w 2 ^i + 



02— Sl 



G 3 (to) + 

1 e*^ 



#2 



+ 



(A4) 



+ 



W 3 ( ^3^-1 e J 2 A '-l \ 
72 ^ 3 _^ 2 I ^3 02 J 



Gs(ti)= a2a3 fe<|«G 1 (t„) + 



+Q3 e % 3 :v G 2 (t ) + 

+ e^ At G 3 (i„) + 

-I-™ ^ r g3(At)- 92 (At) _ g3 (At)- gl (At) 1 

^"^2-*! L <?3-#2 #S-*1 J 

_1_ (y., W 2 / =1 _ 

+ "3 3 -1?2 V *3 #2 

+ (^3+^7)^1 



+ 



With zero flow velocity (a* = ft = 7fe = 0), (A4) reduces to 
the solving formula of the original static model (see PCB98). 
Notice that the solution Gi(ti) for the I s shell includes not 
only the contribution of the contiguous 2 nd shell, but also a 
contribution from the 3 rd shell "scaled" by its passage through 
the 2 nd shell. Similarly, the 3 rd shell is affected not only by the 
2 nd , but also by the 1 st shell though they are not contiguous. 

With analogous procedure, for an arbitrary number N of 
shells the solution is of the kind: 



Gk(t) = *w(At) Giito) + H ^(At) w " 



(A5) 



where: 












F kk (At) = 


e"* At 










Fki(At) = 


7fc 7fc+i- 


H 


-1 J~k(k + \).. 


AAt) 


I > k 


FkiiAt) = 


Qi + l 


■ Ctk-l 


ak^m+i)... 


AAt) 


I < k 


H kk iAt) = 


ffk(Ai) 










H k iiAt) = 


7fc lk+i- 


H 


-lHk(k+X).. 


AAt) 


I > k 


ffw(At) = 


Ctl + l 


■Ctk-l 


akTi.k(k-i) 


..AAt) 


I < k 
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and the quantities T and TC are constructed by means of re- 
cursive formulae: 



F ki (At) = f ki (At) 

•Tkij — 



J-kijm — 



H 



kij 



kijrn 



Fkijm "Phi 



kijmn — ij m — 19„ 



fit H-kijm ^-kijn 

rtkijmri dm—$n 



The coefficients Fki and Hki describe the contribution of the 
generic shell I to the chemical evolution of k. Notice that a 
shell external to k (I > k) can influence k only if all the inflow 
coefficients 7 in between / and k are non-zero, namely if there 
is a continuous inflow from I to k, as expected. The same holds 
for inner shells I < k, whose contribution Fki, Hki is non-zero 
only if none of the intermediate outflow coefficients a is zero. 
The solution ( |A5| ) gets more and more complicated the larger 
the number N of shells considered, since each shell formally 
feels the contribution of all the other shells (as already noticed 
in the above case N = 3). This occurs because ( |AE| ) would be 
the exact analytical solution in the case of a differential system 
with constant coefficients, namely if the A matrix in (^) were 
constant. Then, (A5) would describe the complete evolution of 
any shell k, which over a galaxy's lifetime would indeed pro- 
cess and exchange gas drifting from or to rather distant shells. 
But A(t) is not constant even when the flow pattern a, j3, 7 is 
constant, because r/k{t) a nd th erefore evolve in time due 

to SF; in fact, we apply (A5) only upon short timesteps At, 
within which A(t) can be considered approximately constant. 
If At is short enough with respect to the radial flow veloci- 
ties — as guaranteed by the Courant condition At < vAr, see 
S 3.2 — , we can assume that within At the fc-th shell is af- 



fected just by the flows from the contiguous shells fe+1 and 
fc-1, and not from more distant shells, although all of them 
formally contribute to the solution. In this approximation we 
neglect all higher order terms in a and 7, namely all the terms 
O(0iiOLj) and 0(7;7j), keeping only the "linear" terms of the 
kind akfk(k-i) and Jkfk(k+i)> so the general solution (A.5) re- 
duces in fact to (12(1) . 



Appendix B: Testing the numerical model 

Since discretized numerical solutions for partial differential 
equations containing an "advection term" tend to be affected 
by instability problems (e.g. Press et al. 1986), we tested our 
numerical code with gas flows against suitable analytical cases. 
We report here two representative tests involving pure gas flows 
(no SF) which allow for exact analytical solutions; to these we 
compare the predictions of our numerical code with a SF effi- 
ciency dropped virtually to zero. 

Selecting the timestep. The first reference analytical case is 
that of a purely gaseous disc of infinite radial extent and flat 
profile, where the gas is: 
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Fig. Bl. Numerical models compared to the exact ana- 
lytical solution for a flat accretion profile with inflows. 
Different models correspond to different typical timesteps 
(see legend on top right) and to different shell spacing 
(upper to lower panel). 



2. flowing inward with a constant (in time) and uniform (in 
space) velocity v, starting from the "inflow onset time" T r /. 

Such a system is governed by the same differential equa- 
tion @ describing the adopted bound ary co ndition at the 
outer disc edge in the chemical model (§3.1.2). Eq. dig) ) is a 
linear, first order, partial differential equation, equivalent to 
the system: 



da 
dt 



Ae~ 



(a) 
-° (b) 



(Bl) 



Eq. (Bla) is solved as: 



r =v(t - to) + r 



and substitution into Eq. (Bib) yields: 



1. accreting uniformly with a time-scale r; 



da 1 . _ 

— H a = Ae 

dt t+p 



P = to 

V 



(B2) 
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This linear, first order, ordinary differential equation is solved 
as: 



r(t) =a(to)e 5 + pde + f* e 



Ae--d£ 



(to + p) e~~ - (t + p) e~~ + r ( e~~ 



Finally, replacing back: 
ro = r - v (t - to) , 
we get: 



t + p = - 

v 



y{r,t)= [l-^(t-to)]cr(r-v(t-to),to) + At x 



?(*-*>)) 



e t — e 



+ -re 



(B3) 



Let's now set the initial conditions at to. If radial flows "ac- 
tivate" at a time to = T r f > the surface density distribution 
for t < T r f is determined just by the accretion profile: 



a(r,T, 



rfj 



At 1 



(see Eq. ^) and Eq. 



^B3|) becomes in fact Eq. ([jT]). Here in 
our test case, Eq. (tL4) is the exact analytical description of 
the surface density profile over the whole disc (a part from the 
centre r = 0, which is a singular point). 

As a representative test, let's consider the case r = 3 Gyr, 
T r f — and v = —1 km sec . The relevant analytical solution 
is plot in Fig. Bl for t = tG = 15 Gyr (thick solid line). The 
numerical models used for comparison cover the radial range 
2 to 20 kpc and adopt a flat accretion profile A(r) = A. Their 
outer edge does match exactly with the analytical counterpart, 
since the boundary condition at r = 20 kpc is given by the an- 
alytical expression ( |l4|) itself. But the predictions of numerical 
models at inner radii tend to deviate from the reference density 
profile, and the mismatch is larger: 

1. the larger the typical timestep of the model; 

2. for a fixed timestep, the thinner the shells (compare upper 
to lower panel). 

For a typical shell width of 1 kpc, for instance, a good match 
is obtained with model timesteps of 2 x 10 -2 Gyr, while model 
shells of 0.5 kpc an acceptable profile is obtained only with 
timesteps of 2 x 10 -3 Gyr. Therefore, a reliable representa- 
tion of radial gas flows is obtained only with a suitably small 
timestep; how small, is related to the width of the shells, 
namely to the resolution of the grid spacing. 

Selecting the grid spacing. Since our models are to simulate 
a disc with an exponential profile, as a second test we consider 
a gaseous disc with uniform and constant infall time-scale and 
inflow velocity, analogous to the previous case, but with an 
exponential profile. Namely, in the representative differential 
equation (113) the accretion profile A(r) declines exponentially 
outward: 




Fig. B2. Exact analytical solution for an exponential ac- 
cretion profile with inflows compared to a numerical model 
with 40 shells equally spaced in logarithmic scale. 

Eq. (N 2h can then be written: 



A{ro) e r d e p ie 



da da 

This is another linear, first order, partial differential equation 
of the same kind as (p^[), and can be solved with analogous 
procedure into: 

a{r,t) = [l-£.(t-t )]ff(r-v(t-t ),to) + 

[(l-H*-*o))e-^ (t - to) + (B4 ) 



-e~~ + 



Notice that Eq. (B3) for a flat profile is recovered from (B4) 



for rd — ► oo. If radial inflows set in at time to = T r f > 0, 
from 



a{r,T rf ) =A{r Q ) e ^ r (1 - e - 
we get: 

a(r, t) — A(r Q ) e r a ' |r ^1 — e~~ 

-K-t— ( i + )T rf t 



+ 



X t|1 - e 



(B5) 



- (t-Trf)x 



A(r) = ^(ro) 



We compared this second analytical case to our numeri- 
cal model, where an exponential accretion profile was adopted 
between 2 and 20 kpc, and at the outer edge the boundary con- 
dition (O) was replaced by ( B5). Our tests showed that in this 
case the analytical profile is better reproduced the larger the 
number of shells, i.e. the finer the grid spacing. A good match is 
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obtained especially when model shells are equally spaced in the 
logarithmic, rather than linear, scale; namely when the shells 
are chosen, in this case of an exponential accretion profile, so as 
to contain roughly the same mass, rather than cover the same 
radial width. Fig. B2 shows in fact the analytical solution (B5) 
(for Td = 4 kpc, r = 3, T r f = and v = — 1) together with a 
corresponding numerical model with 40 shells logarithmically 
spaced (from 0.1 kpc wide for the inner ones to ~1 kpc wide 
for the outer ones). 

With such a grid spacing, our tests with the reference flat 
indicate 10 -4 Gyr as a suitable timestep 



Bl 



profile of Fig. 

to obtain stable solutions for velocities up to the order of 
1 km sec - . 

In the light of all these tests, in our chemical models we 
adopted a grid spacing of 35 shells from 2.5 to 20 kpc, equally 
space d in logarithmic scale and a typical timestep of 10 -4 Gyr 
(see §3.2). Of course, the suitable timestep depends also on the 
velocity field: flows with higher velocities require shorter inte- 
gration timesteps. Whenever we need to consider much larger 
speeds than 1 km sec -1 , as might be the case for the strong 
flows induced by the Bar, we reduce the timestep in proportion. 
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